Direct observation of a Berezinskii-Kosterlitz-Thouless superfluid in an atomic gas 
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Understanding the emerging mechanisms of superfluidity has been a central theme in many-body 
physics. In particular, the superfluid state in a two-dimensional (2D) system is intriguing because 
large thermal fluctuations prohibit the formation of long-range order [TJ [2] and consequently the 
picture of Bose-Einstein condensation is not applicable to the phase transition. The Berezinskii- 
Kosterlitz-Thouless (BKT) theory [3j [4] provides a microscopic mechanism for the superfluid phase 
transition, where vortices with opposite circulation are paired below the critical temperature, es- 
tablishing quasi- long-range phase coherence. While this mechanism has been tested experimen- 
tally [SHS], there has been no direct observation of thermal vortex pairs in a 2D superfluid. Here 
we report the observation of a BKT superfluid in a trapped quasi-2D atomic gas by measuring the 
spatial distribution of thermally activated vortices and revealing their pair correlations. The vortex 
population concentrates in the low-density superfluid region of the trapped sample and decreases 
gradually as the temperature is lowered, showing a crossover from a BKT phase to a Bose-Einstein 
condensate (BEC) with no thermal vortices. These observations clarify the nature of the superfluid 
state of a trapped 2D Bose gas. 



The BKT mechanism represents a topological phase 
transition involving no spontaneous symmetry-breaking. 
The pairing of vortices is responsible for the emergence 
of a topological order for the low-temperature superfluid 
phase. Since vortex- ant i vortex pairs carry a zero net 
phase slip, the BKT superfluid can be described to be 
topologically identical to a state with a spatially uni- 
form phase like a BEC on the large length scale com- 
pared to the vortex pair size, although it contains long- 
wavelength phase fluctuations due to thermal phonon ex- 
citations. Above the critical temperature, the topologi- 
cal character of the system disappears as free vortices, 
i.e. point defects, proliferate via unbinding of the vortex 
pairs, which is reflected in the decay behavior of phase 
coherence changing from algebraic to exponential. The 
BKT mechanism has been supported by experimental ob- 
servations in many 2D systems [5H9]. However, a direct 
observation of thermally activated vortex pairs has been 
elusive. 

Ultracold atomic gases in 2D geometry present a clean 
and well-controlled system for studying BKT physics. In 
particular, individual vortices can be detected in this sys- 
tem, providing interesting opportunities to directly test 
the microscopic nature of the BKT transition. Previ- 
ous experimental studies of phase coherence [8] [10] and 
thermodynamic properties [TTHT3] indicated that a BKT- 
type transition occurs in a quasi-2D Bose gas trapped 
in a harmonic potential. Recently, superfluid behavior 
was demonstrated by measuring a critical velocity for 
friction- less motion of an obstacle [14]. In the present 
study, we measure the in-plane spatial distribution of 
thermally activated vortices in a quasi-2D Bose gas and 
reveal the vortex pairing from their spatial correlations. 
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In addition, we observe that a vortex-free region emerges 
in the center of a trapped sample below a certain temper- 
ature, demonstrating crossover behavior of the superfluid 
state from a vortex-pair BKT phase to a vortex-free BEC 
phase [15j [16] . Our results clarify the nature of the su- 
perfluid state in an interacting 2D Bose gas trapped in a 
harmonic potential. 

Our experiments are carried out with a quasi-2D Bose 
gas of 23 Na atoms in a pancake-shaped optical dipole 
trap [17 . A highly saturated thermal gas of 23 Na atoms 
in the \F = l,m^ = —1) state is first generated in an 
optically plugged magnetic quadrupole trap, and trans- 
ferred into an optical dipole trap that is formed by focus- 
ing a 1064 nm single laser beam. The final evaporative 
cooling is achieved by lowering the trap depth for 9 s. 
The trapping frequencies of the harmonic potential are 
(uj x ,ujy : uj z ) = 2tt x (3.0,3.9,370) Hz, where the z-axis is 
along the direction of gravity. The dimensionless interac- 
tion strength g = \^S7ra/l z c± 0.013, where a is the three- 
dimensional (3D) scattering length and l z = ^/h/muj z 
is the axial harmonic oscillator length (m is the atomic 
mass). For a typical sample of N = 1.3 x 10 6 atoms, 
the zero-temperature chemical potential is estimated to 
be fx ~ 0.7Hoj z , where H is the Planck constant h divided 
by 27r, and thus, our sample constitutes a very weakly 
interacting quasi-2D Bose gas. 

In this harmonic trap, the Bose-Einstein condensa- 
tion temperature for a quasi-2D ideal gas of N atoms 
is Tq ~ 76 nK (To is calculated from the relation N = 

h 2 u uj ^2^=0 Li2(e Z/fcT °), considering thermal popula- 
tion in the tight direction, where Li 2 (z) is the polylog- 
arithm of order 2) and the BKT critical temperature is 
estimated in the mean- field theroy to be T c /Tq ~ 0.99 
for our small g = 0.013 [18]. Below the superfluid phase 
transition, the sample shows a bimodal density distribu- 
tion after time-of-flight expansion (Fig. la) [19 . We refer 
to the center part as the coherent part of the sample and 
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FIG. 1: Observation of thermally activated vortices, (a) Opti- 
cal density image of a quasi- 2D Bose gas, taken after a 17-ms 
time-of-night expansion. Density ripples develop during the 
expansion due to phase fluctuations in the sample, (b) Im- 
age with applying a radial compression before the expansion. 
The radial compression is achieved by superposing a magnetic 
potential onto the sample in an optical trap. Quantized vor- 
tices are observed with density-depleted cores, (c) Positions 
of vortices with respect to the sample center recorded for 20 
experimental realizations. The average number of vortices ob- 
served in a sample is N v ^17. (d) Radial profile of the vortex 
density n v (r). The red solid line displays the atom density 
profile n(r) of the sample, showing a bimodal distribution. 
The center part is refered to as the coheret part of the sample 
and its boundary is indicated by an arrow in (d) and a red 
solid line in (c). The data are acquired for N = 1.6(2) x 10 6 
atoms at T = 50(6) nK. 



determine the sample temperature T from the gaussian 
width of the outer low-density thermal wing. 

The conventional method for detecting quantized vor- 
tices is observing density-depleted vortex cores after re- 
leasing the trap [20]. Because the vortex core size £ ex 
n -1 / 2 oc i? 3 / 2 , where n is the atom density and R is 
the sample size, the vortex core expands generally faster 
than the sample in a 3D case [21 , facilitating its detec- 
tion. For a 2D Bose gas, however, this method is not 
adequate because of the initial fast expansion along the 
tight confining direction which reduces atom interaction 
effects so rapidly that the vortex core structure is not 
maintained during the subsequent evolution. Further- 
more, phase fluctuations present in the 2D sample evolve 
into density modulations [17 , making the detection of 
individual vortices more difficult. 

To improve the visibility of the vortex cores, we apply a 
radial compression to the sample before expansion. The 
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FIG. 2: Pair correlations of vortices, (a) Two- vortex spatial 
correlation function g(s), obtained from 20 images for the ex- 
perimental condition in Fig. 1. The red solid line represents 
the correlation function go(s) for an uncorrelated distribution 
with the same vortex density profile n v (r) in Fig. lc. The in- 
set displays Fig. lb with guide lines for pairs of vortices whose 
separation is less than 30 /mi. (b) The difference g(s) — go(s) 
shows oscillatory behavior with a noticeable enhancement at 
s ~ 35 /im. 



compression transforms the 2D sample into an oblate 3D 
sample with /jl/Hlj z > 1, restoring the favorable condition 
for the vortex detection [21] . Furthermore, in the com- 
pressed sample, thermal relaxation of phonon excitations 
will proceed because of the change in the dimensionality 
of the sample. Vortices may also decay, but the effective 
2D geometry in terms of vortex relaxation can be pre- 
served for a moderate compression, suppressing Kelvin 
mode excitations on a vortex line [22] . 

In our experiment, the radial compression is achieved 
by superposing a quadrupole magnetic potential onto the 
sample. The axial axis of the magnetic quadrupole field is 
aligned with that of the sample. We control the field gra- 
dient and the position of the zero-field point to increase 
the trapping frequencies at the center of the hybrid trap 
almost linearly to uo x ^ y = 2tt x 39 Hz for 0.4 s, and let 
the sample relax for 0.2 s that corresponds to about 10 
collision times. There was no significant collective oscilla- 
tion of the sample, and the atom number fraction for the 
cohrent part was reduced by about 7% in this process. 
The chemical potential of the compressed sample is esti- 
mated to be ji/hw z > 5, reaching a 3D regime. This was 
confirmed by making a condensate in the tight trap and 
observing no phase fluctuations [17]. To magnify vortex 
cores, we employ a two-step expansion: First, we sud- 
denly turn off only the magnetic potential to make the 
sample expand radially for 6 ms in the optical trap, and 
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FIG. 3: Compression time dependence of the vortex distribution, (a) Evolution of the scaled chemical potential p/huj z in the 
radial compression process, (b) Vortex number N v and (c) nearest vortex distance d m versus compression time r. (d, e, f) 
Radial profiles of the vortex density for various compression times. The vortex density profile for r = 0.4 s is displayed in the 
background and the red solid lines show the atom density profiles. These data are acquired for N — 1.8(2) x 10 6 atoms at 
T — 46(6) nK. Error bars indicate standard deviation. 



then switch off the optical trap for further magnification. 
It is not clear why the in-trap expansion could enhance 
the contrast of vortex cores, which is not an exact 2D 
expansion. Finally, we take an absorption image along 
the axial direction after additional 9 ms of time-of-flight 
expansion to measure the in-plane density distribution of 
the sample. 

Using this procedure, we observe thermally excited 
vortices with clear density-depleted cores (Fig. lb). This 
indicates the existence of vortex-antivortex pair excita- 
tions because the direct generation of single vortices is 
forbidden. The core size is measured to be £ v ~ 8 /im, 
allowing by-hand locating of the vortices that gives the 
vortex distribution p v (r) = E^5(r — f$), where N v is 
the number of voritces and ?i is the position of the i- 
th vortex with respect to the center of the sample. At 
T « 50 nK, the average vortex number is about 17 and 
vortices mainly appear in the outer region of the coher- 
ent part. There is no azimuthal dependence of the vortex 
distribution (Fig. lc) and the radial profile of the vortex 
density n v {r) is obtained by azimuthally averaging the 
averaged p v (r) for the same experiments (Fig. Id). The 
vortex region (n v > 0) reaches the boundary of the co- 
herent part, implying the existence of a vortex-driven 
phase transition. Density modulations suggestive of vor- 
tex cores are often observed at the boundary, but not 
included in the vortex counting. 

One notable feature in the vortex distribution is fre- 
quent appearance of a pair of vortices that are closely 
located to each other but well separated from the others 
(Fig. 2 inset). We measure the distance from every vor- 
tex to its nearest vortex and find its average value to be 
dm ~ 4£ v . A more quantitative study of the pair correla- 
tions is performed with the two- vortex spatial correlation 



function: 

«'>- ^.w-i) g"- r «> (1> 

l>j 

where = \fi — fj\. This function displays the prob- 
ability of finding a vortex at distance s from another 
vortex, containing rich information about vortex inter- 
actions. We determine g(s) as the average of the pair 
correlation functions individually obtained from images 
for the same experiments (Fig. 2a). In order to see 
vortex interaction effects, we compare g(s) to the pair 
correlation function for an uncorrelated vortex distri- 
bution having the same vortex density profile n v (r): 
go(s) = j drdO rn v (r)n v (r f )/[j dr27rrn v (r)] 2 where r' = 
VV 2 + 2rs cos + s 2 . The difference g(s) — go(s) shows a 
noticeable enhancement around s ~ dm and small oscil- 
latory behavior at large s, clearly indicating that there 
are attractive interactions between closely located vor- 
tices. The complete suppression at s < 3£ v is attributed 
to the annhilations of tightly bound vortex pairs during 
the compression and expansion process. When we slowly 
compress the sample, the vortex number N v rapidly re- 

— 1/2 

duces over 1 s, but d m increases more slowly than N v 
(Fig. 3b and c). 

The pairing feature in g(s) becomes more pronounced 
at lower temperatures (Fig. 4i-k), which we ascribe to 
the temperature dependence of the characteristic size of 
vortex pairs. The thermal excitation probability p of 
a vortex-antivortex pair is estimated from the associ- 
ated free energy F = E — TS as p oc e~ F l kT where 
k is the Boltzmann constant. The excitation energy 
E of a vortex pair of size d is ~ ^^n s \n(d/^) where 
n s is the superfluid density [15 j. In a superfluid re- 
gion of size R, the number of distinguishable microstates 
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FIG. 4: Thermal excitations of vortices in a trapped 2D Bose gas. The left column (a-d) displays examples of vortex images for 
various temperatures. The experimental parameters and the votex number are (N,T,N V ) = (1.5(2) x 10 6 , 67(12) nK, 13(4)) 
for (a), (1.4(3) x 10 6 ,56(6) nK, 16(4)) for (b), (1.3(1) x 10 6 ,48(8) nK, 15(5)) for (c), and (1.0(3) x 10 6 ,23(4) nK, 2(4)) for (d). 
The middle column (e-h) and the right column (i-k) show the radial profiles of the vortex density n v (r) and the pair correlation 
functions g(s), respectively, for the experimental conditions in (a-d). The data are obtained from at least 12 images for each 
condition. g(s) for the lowest temperature (d) is not reliably determined due to the small vortex number. The red solid lines 
in (e-h) and (i-k) represents n(r) and go(s), respectively. 



for the vortex pair is ~ f^^f > where the first and the 
second terms correspond respectively to the position of 
a single vortex in the region and the relative orien- 
tation of two vortices with separation d, which gives 
the entropy S = k ln(7rdR 2 /£ 3 ). In this estimation, 
p oc ^(d/0" (nsA2_1) where A = h/V^TrmkT is the ther- 
mal wavelength. Therefore, tightly bound pairs are rela- 
tively more populated at lower temperatures, leading to 
more enhancement at smaller s in the pair correlation 
function. These obsevations provide conclusive evidence 
for the pairing of vortices in the superfluid phase in a 
quasi-2D Bose gas. 

Figure 4e-h display the vortex density profiles n v (r) for 
various temperatures. Because the thermal excitation 
probability is exponentially suppressed by n s A 2 , n v (r) 
reflects the spatial distribution of the superfluid density 
n s (r). At high temperature near the critical point, vortex 
excitations prevail over the whole coherent part, and as 
temperature is lowered, a vortex-free region starts emerg- 
ing in the center of the sample, pushing vortices out- 
wards. At T « 30 nK, only a few vortices appear near 
the boundary of the coherent part. The emergence of 
a vortex-free region indicates the change in the micro- 



scopic nature of the superfluid phase. This represents a 
crossover from a BKT phase to a BEC phase, continously 
suppressing the vortex-pair excitations. In a finite-size 
2D system, forming a BEC is expected at finite tempera- 
ture when the coherence length of the superfluid becomes 
comparable to the spatial size of the system [I5j [16] . 



The BKT theory for a homogeneous system predicts 
spontaneous unbinding of vortex pairs and a sudden drop 
of the superfluid density at the critical point [23] . These 
would be revealed in our experiments as a rapid change in 
the pair correlations or/and sharp behavior of the outer 
boundary of the vortex region. However, these are not 
observed in our experiments. This might be due to the 
limited sensitivity of our detection method in the crit- 
ical region and possibly, drifting of vortices out of the 
coherent part during the compression process. On the 
other hand, the distinction between vortex pairs and free 
vorices is not clear in a finite system and the behavior 
of vortices might be expected to be continuous over the 
critical point [8]. 

We summarize our results in Figure 5 with a phase di- 
agram of a 2D Bose gas in a harmonic potential. The in- 
situ critical radius R c for the superfluid-to-normal phase 
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FIG. 5: BKT-BEC crossover. The phase diagram of a 2D 
Bose gas trapped in a harmonic potential is depicted in the 
plane of r / Rtf and T/Tq. Rtf is the 2D Thomas- Fermi ra- 
dius, and To is the BEC temperature for a quasi- 2D ideal 
gas in the harmonic potential. A BKT superfluid phase con- 
taining vortex-pair excitations emerges below a critical tem- 
perature. The in-situ critical radius R c for the superfluid- 
to-normal phase transition are indicated by open circles, 
estimated from the radius of the coherent part, including 
the time-of- flight expansion factor. The red dashed line is 
v/Rtf = yl — (T/Tb) 4 - 2 as a guide for the critical line. As 
the temperature is decreased, thermal excitations of vortex 
pairs are suppressed, preferentially in the center region of high 
superfluid density, leading to the emergence of a BEC phase at 
low temperature. Red solid circles represent the inner radii 
R v of the vortex region (n v > 0), determined with the as- 
sumption that R v /R c is preserved during the compression and 
expansion process. Error bars indicate standard deviation. 

transition is determined from the radius of the coherent 
part in the x-direction measured in images taken without 
the compression process as Fig. la, including the expan- 
sion factor for 17-ms time-of- flight that was measured to 
be 1.03(1) from a fit of R c (l — at 2 ) to the radii of the 



coherent part for various time-of- flight t = 10 ~ 25 ms. 
The 2D Thomas-Fermi radius R TF = f 
the x-direction. The vortex region is marked as a BKT 
phase, where the inner radius R v is determined from a 
hyperbolic-tangent fit to the inner region of n v (r) with 
a threshold value of 6 x 10 -5 fim~ 2 . Here, we assume 
that the relative radial position R v /R c is preserved in 
the radial compression and expansion process. Because 
a vortex pair carries a linear momentum of h/d, the vor- 
tex region might be broadened by diffusion during the 
compression. However, we observed that the inner radius 
R v of the vortex region does not change significantly for 
longer compression times (Fig. 3d-f), probably because 
it is energetically unfavorable for vortices to intrude into 
the center region with high superfluid density. In re- 
cent theoretical studies [24, 25 , the critical temperature 
for the emergence of a BEC phase was calculated as a 
point when the thermal excitation probability of a vortex 
pair becomes of order unity. Our results show qualitative 
agreement with the predictions, but their direct compar- 
ison is limited because the vortex detection efficiency of 
our imaging method is not determined. 

In conclusion, we have observed thermally activated 
vortex pairs in the superfluid phase of a 2D Bose gas 
trapped in a harmonic potential. This provides the clear 
confirmation of BKT superfluidity in the systems. The 
experimental study of 2D superfluidity has been recently 
extended to ultracold atom systems including disorder 
potentials [26] and Fermi gases [27] , and new research 
directions with spin degrees of freedom [28] and in a dy- 
namic transition regime [29] have been theoretically pro- 
posed. The vortex detection method developed in this 
work will be an important tool to probe the microscopic 
properties of these systems. Integrated with an intefero- 
metric technique [10] [30], this method might be upgraded 
to be sensitive to the sign of a vortex. 
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